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We present the first measurement of T(4S') decays to 7r"'"7r T{1S) and ■k'^-k T{2S) based on 
a sample of 230x10® T{AS) mesons collected with the BABAR detector. We measure the product 
branching fractions B(r(4S') -^ 7r+7r"r(lS))xS(r(lS') -^ m^M") = (2.23±0.25^tat±0.27sas)x lO"*^ 
andS(r(45) -* ■K+Ti-T{2S))xB{Y{2S) -> m^m") = (1.69±0.26stat±0.20sj,s)xl0-®, from which we 



derive the partial widths r(T(45') 
(2.7 ±0.8) keV. 

PACS numbers: 14.40.Gx,13.25.Gv 



-T{IS)) = (1.8±0.4) keV and r(r(45') 



-T{2S)) 



The T{AS) meson is known to decay predominantly to 
BB, with small, but as yet unobserved, decays to other 
bottomonium states or to light hadrons. Partial widths 
for hadronic transitions in heavy quarkonia have been ex- 
tensively studied both experimentally and theoretically 
over the past decades []|. In particular, the values of 
the partial widths for dipion transitions between vector 
states ip{2S) -^ tt+tt-J/iP and T{mS) -^ Tr+Tr~r (nS), 
where the principal quantum number m > n, can be re- 
lated to the radial wave function within the framework 
of the QCD multipole expansion [2|. This picture may 
be significantly altered by mixing and coupled channel 
effects |3| when states are close to the threshold for open 
charm or bottom production. Hence these states are the 
ideal laboratory to investigate these effects. Exclusive 
non-D_D decays of the ^"(3770) (believed to be predom- 
inantly '^Di) have recently been observed |j, laldi but 
only upper limits have been published for exclusive non- 
SB decays of the r{4S) 0. 

We search for the decays T(4S') -^ 7r+7r~T(n5'), where 
n = 1, 2 0, using a sample of 230 x 10^ T{4S) events 
corresponding to an integrated luminosity of 211 fb~^ ac- 
quired near the peak of the T{AS) resonance with the 
PEP-II asymmetric-energy e~^e~ storage rings at SLAC. 
An additional 22 fb~^ sample collected approximately 
40 MeV below the resonance is used as a control sam- 
ple. 

The BABAR detector is described in detail elsewhere |9| ; 
here we summarize only the features relevant to this anal- 
ysis: charged-particle momenta are measured in a track- 
ing system consisting of a five-layer double-sided silicon 
vertex tracker (SVT) and a 40-layer central drift cham- 
ber (DCH), both situated in a 1.5-T axial magnetic field. 
Charged-particle identification is based on the dE/dx 
measured in the SVT and DCH, and on a measurement of 
the photons produced in the synthetic fused-silica bars of 
the ring-imaging Cherenkov detector (DIRC). A CsI(Tl) 
electromagnetic calorimeter (EMC) is used to detect and 
identify photons and electrons, while muons are identified 
in the instrumented flux return of the magnet (IFR). 

An T{mS) —> TT^n^T{nS) transition, denoted by 
mS -^ nS, is detected by reconstructing the T{nS) me- 
son via its leptonic decay to fJ.^ n~ ■ The sensitivity to 
45 — > nS transitions is much smaller in the 7r+7r~e^e^ 
final state due to the presence of larger backgrounds, and 



to a trigger-level inefficiency introduced by the prescaling 
of Bhabha scattering events. Data collected at a nom- 
inal center-of-mass energy y/s near 10.58 GeV include 
35 — > nS {n — 1, 2) and 25 — > 15 events from initial 
state radiation (ISR) production that are used as control 
samples. The signature for mS —f nS transition events, 
where the nS decays to muons, is a fJ-^ fJ-~ invariant mass, 
Af^^, that is compatible with the known mass [lO| of 
the Y{nS) resonance, M{nS), and an invariant mass dif- 
ference AAf = A/^TrpM ^ -^'^MM ^^^^ is compatible with 
M{mS) — M{nS). The r.m.s. values of the reconstructed 
AM and M^^ distributions are, respectively, « 7 MeV/c^ 
and ~ 75 MeV/c^. The center-of-mass momentum p*^^j^ 
should be compatible with for 45 — > nS candidates, 
or with (s — M'^{mS)) /{2y/s) for mS — > nS candidates 
from ISR. 

Simulated Monte Carlo (MC) events are generated us- 
ing the EvtGen package [ll|. The angular distribution of 
generated dilepton decays incorporates the T{nS) polar- 
ization, while dipion transitions are generated according 
to phase space. These events are passed through a detec- 
tor simulation based on GEANT4 [l2|, and analyzed in 
the same manner as data. The events in the data sam- 
ple whose values of AM and M^^ are within 60 MeV/c^ 
and 300 MeV/c^ , respectively, of the values expected for 
any known mS — > nS transition were not examined until 
the event selection criteria were finalized. Events outside 
these regions were used to understand the background. 

We select events having at least 4 charged tracks with 
a polar angle 9 within the fiducial volume of the tracking 
system (0.41< 9 <2.54 rad). Each muon candidate is re- 
quired to have a center-of-mass momentum greater than 
4 GeV/c, and to be compatible with the muon hypothesis 
based on the energy deposited in the EMC and the hit 
pattern in the IFR along the track trajectory. A dipion 
candidate is formed from a pair of oppositely charged 
tracks. The two pion candidates are each required to 
have a transverse momentum greater than 100 MeV/c. 
The dimuon and the dipion are constrained to a com- 
mon vertex, and the vertex fit is required to have a x^ 
probability larger than 10""^. 

A large fraction of the background is due to fj,'^fj.~j 
events where a photon converts in the detector mate- 
rial. To reduce this background we apply an "elec- 
tron veto", rejecting events where any of the follow- 
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limit the regions where AM and M^^ are within ±60 MeV/c'^ 
and ±300 MeV/c^, respectively, of the values expected for an 
mS —> nS transition. The remaining region is used to model 
background. The text discusses the features seen in the data. 



ing is true: either of the two pion candidates is pos- 
itively identified as an electron; the e"*"e~ invariant 
mass of the two charged tracks associated with the 
pion candidates satisfies M^e < 100 MeV/c^; or the di- 
pion opening angle satisfies cos^tt+tt- > 0.95. The 
distribution of AM vs M^^ for the final sample is 
shown in Fig. ^ The clusters of events in the boxes 
centered at (AM,M^^) = (1.120, 9.460) GeV/c^ and 
(0.558, 10.023) GeV/c^ constitute, respectively, the first 
observation of 45 — » IS* and of 45 -^ 25 transitions. We 
also observe signals for 25 -^ 15, 35 -^ 25, and 35 -^ 
15 from ISR at (AM,M^^) = (0.563,9.460) GeV/c^, 
(0.332, 10.023) GeV/c2, and (0.895, 9.460) GeV/c^ respec- 
tively. The diagonal band is predominantly due 
to fifi'j events, while the cluster at (AM, M^^) = 
(0.332,9.460) GeV/c^ is due to r(35) -^ ti+tt-T{2S) de- 
cays, where r(25) ^ r(15) X. 

The number of signal events Nsig is extracted by an 
unbinned extended maximum likelihood fit to the AAf 
distribution for events with p*^a,nd < 200 McV/c and 
\M^^ - M(15)| < 200 MeV/c^ for the 45 -^ 15 mode 
or |M^^ - M(25)| < 150 MeV/c^ for the 45 ^ 25 mode 
(Fig. 12) ■ In each case, the background is parametrized 
as a linear function, and the signal as the convolution 
of a Gaussian with standard deviation a and a Cauchy 
function with width F, which is found to adequately de- 
scribe the non-Gaussian tails of the AM distribution. 
The values for a and F are, for each mode, fixed to 
the values determined from a fit to a MC signal sam- 
ple subjected to the detector simulation and reconstruc- 
tion algorithms. We verify that the experimental AM 
resolution is well described by the MC for 25 ^- 15 
and 35 — > nS {n = 1, 2) ISR samples. The values of 
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FIG. 2: The AM distribution for events with |M^+^- - 
M(1S')| < 200MeV/c2 (left) and |M^+f- - Af(25')| < 
150MeV/c^ (right). The solid hues show the best fit to the 
data. Dashed lines show the background contribution. The 
two upper plots are for tv'^ tt' /i^ fi~ candidates and the two 
lower plots for 7r"'"7r~e^e~ candidates. 



AM returned by the fit, 1.1185 ± 0.0009 GeV/c^ and 
0.5571 ± 0.0010 GeV/c^, where the errors are statisti- 
cal only, are in excellent agreement with the world av- 
erages M(45) - M(15) = 1.1197 ± 0.0035 GeV/c^ and 
M(45) - M(25) = 0.5567 ± 0.0035 GeV/c^ |i(j]. These 
values cannot be interpeted as a new measurement of 
the T(45) mass, since data come from the peak of the 
resonance and not from a scan of it. The cuts described 
above are also applied to 7r"'"7r~e"'"e^ candidates, with the 
additional requirement on the polar angle of the electron, 
0{e~) > 0.75 radians, to reject Bhabha events. The fits 
to the electron samples are also shown in Fig. [2 and give 
yields and AM values consistent with expectations based 
on the fits to the muon samples. 

The significance, estimated from the likelihood ratio 
na ~ ^y2\og[£{Nsig) / £{0)] between a fit that includes a 
signal function and a fit with only a background hypoth- 
esis, is IO.Oct for the 45 -^ 15 and 7.3ct for the 45 -^ 25 
in the 7r+7r^/^+/i^ final states. The significance of the 
signals in the 7r+7r^e+e~ final states is 3.6cr and 2.5(T for 
45 -^ 15 and 45 -+ 25, respectively. 

The event selection efficiency esei is determined using 
the MC samples. The largest source of systematic un- 
certainty (10%) is due to the unknown distribution of 
the dipion invariant mass in the T(45) — » 7r+7r^T(ri5) 
transition, and is estimated by comparing the acceptance 
for a phase space distribution to that obtained using the 



TABLE I: Number of signal events, significance, efficiency 
and measured values of the products of branching ratios for 
the 45* — > nS transitions. The error on the efficiency is 
obtained adding in quadrature the systematic uncertainties. 
The errors on the product branching fractions are statistical 
and systematic respectively 



Transition Nsig significance 



Ssel S4S-»nS X BnS^^l.^i 

(%) (10"^) 



45- 
AS 



IS 
25 



f67±f9 
97±f5 



fO.Ocr 

7.3(T 



32.5±3.9 
24.9±3.0 



2.23±0.25±0.27 
f.69±0.26±0.20 



QCD multipole model '2|. The second largest source 
of systematic uncertainty is due to uncertainty in the 
track reconstruction efficiency, which is 1.3% per track, 
resulting in a 5.2% uncertainty in tsei- The systematic 
uncertainties associated with the event selection (4.3%) 
and muon identification (1.4%) criteria are estimated by 
comparing the efficiency of each selection criterion deter- 
mined from MC samples to the corresponding efficiency 
measured with the ISR control samples. We have also 
considered the systematic uncertainties due to the choice 
of signal and background paranietrizations by using dif- 
ferent functions or different parameters, and the system- 
atic uncertainties due to the choice of the fit range. The 
contributions from these sources are negligible in com- 
parison to the previously mentioned sources. 

The product branching fraction (TableQj is determined 
from the n'^ ■k~ fj,'^ fi~ sample using: 



B{r{AS)-^Tr+Tr-T{nS)) 
B (T{nS) -^ fi+fi- 



^.. 



esei iV(45) ' 



(1) 



where 7V(4S') = (230.0 ± 2.5) x 10^ is the total number 
of T{4S) mesons produced. 

The event yields observed for 3S -^ nS and 25* — + IS 
are compatible with PDG-averaged values of the ISR 
cross section and branching fractions for those reso- 
nances. The number of signal events observed in the 
TT+TT^e+e^ final state is compatible with the branch- 
ing fractions we measure in the 7r+7r^/i+/i^ sample. No 
45 -^ nS signal is observed for ir^ tt~ jX^ ^~ or 7r"'"7r~e"'"e~ 
final states in the data collected at center of mass energies 
40 MeV below the T{AS) resonance. 

The dipion invariant mass distribution, M^+^- 
(Fig. EJ), is determined by fitting the AM distribution 
in equal intervals of M^+^- , and dividing the number 
of signal events in each interval by the corresponding 
selection efficiency. The measured distribution for the 
AS — > IS transition has a shape similar to the prediction 
of the Kuang-Yan model ^ . This model provides a good 
description of the observed distributions for 2S — » IS, 
35' -^ 25, and also ■0(25) -> tt+tt" J/-0, but fails to de- 
scribe the 35 -^ 15 distribution. Our measured distribu- 
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FIG. 3: The efficiency-corrected M^+^- distribution for 
AS -^ IS transition (left) and 4S — > 25' transition (right). 
The solid line shows the distribution predicted in Ref . ,2] . The 
dotted histogram shows the selection efficiency in each bin. 
The experimental resolution in M^+^- is less than 5 MeV/c^, 
much smaller than the bin size. 



tion for the 45 — > 25 transition has a marked enhance- 
ment at low Myr+TT- that is incompatible with this model. 
The 45 — > nS branching ratios and partial widths 
can be derived using the world average values for 
B{T{nS) —> u'^u~) llOl and a recent BABAR measure- 
ment of r(T(45)) [llpWe obtain 



B (T(45) ^ 7r+7r-r(15)) = (0.90 ± 0.15) x 10^^^ 

B (T(45) ^ 7r+7r-r(25)) = (1.29 ± 0.32) x 10-^ 

r(T(45) ^ 7r+7r-r(15)) = (1.8 ± 0.4) keV, 



and 



r(T(45) -> 7r+7r-r(25)) = (2.7 ±0.8) keV. 



We add in quadrature the statistical and systematic un- 
certainties on the derived quantities. With the most 
recent CLEO measurement of 6 (T(25) -^ fi+fi^) Hi, 
we obtain smaller values: K (T(45) — > 7r+7r^T(25)) = 
(0.83 ± 0.16) X 10-4 and r(T(45) -> 7r+7r-T(25)) = 
(1.7 ±0.5) keV. 

The branching fractions are compatible with previ- 
ous upper limits on these decays ^3]. The T(45) partial 
widths are within the raiige spanned by other dipion tran- 
sitions in the bb system "l^: r(T(25) -^ 7r+7r-T(15)) = 
(8.1 ± 2.1)keV; r(r(35) ^ 7r+7r-r(15)) = (1.2 ± 
0.2) keV; r(r(35) -^ 7r+7r-T(25)) = (0.6 ± 0.2) keV. 

In conclusion, we measure 



B (r(45) 



-r(15)) X B (T(15) -^ ^i+fi-) = 
(2.23 ±0.25 ±0.27) x 10"*^ 



and 



B (T(45) ^ 7r+7r-T(25)) x B (r(25) ^ n+ fi') = 
(1.69 ±0.26 ±0.20) X 10"^ 



The dipion invariant mass distribution is measured for 
r(4S) -^ Tr+TT-r{lS) and T{AS) -^ Tr+Tr-T{2S) transi- 
tions; the latter is found to be incompatible with predic- 
tions from QCD multipole expansions. 
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